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<N ■ Abstract 

We revisit the theoretical predictions for anomalous radiative decays of pseu- 
doscalar and vector mesons. Our analysis is performed in the framework of 
the Nambu-Jona-Lasinio model, introducing adequate parameters to account 
^ I for the breakdown of chiral symmetry. The results are comparable with those 

obtained in previous approaches. 
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O : I. INTRODUCTION 

The analysis of vector and axial-vector meson physics represents a nontrivial task from 
P_i' the theoretical point of view. This is basically due to the characteristic energy scales at 
p • which these particles manifest, namely an intermediate range between low energy hadron 
O , physics and the high energy region where QCD can be treated perturvatively. In order to 
deal with the subject, one possible approach is to start from an effective meson Lagrangian; 
another possible way is to derive the corresponding effective interactions from QCD-inspired 
^ ■ fermionic schemes. 

■ In the last two decades, the development of chiral effective models ^ has provided a 

profitable framework to analyze various phenomena related to the physics of vector and 
axial- vector mesons . These models have been built taking into account the symmetries 
of QGD, together with the so-called axial anomaly introduced through the Wess-Zumino- 
Witten (WZW) ||^ effective action. As an alternative approach, the analysis of spin one 
meson physics can be performed by starting with the Lagrangian proposed by Nambu and 
Jona-Lasinio (NJL) ^j. This Lagrangian is based in 4-fermion interactions, showing the 
chiral symmetry of QCD, and leads to an effective theory for scalar, pseudoscalar, vector 
and axial-vector mesons after proper bosonization. The effective meson Lagrangian can be 
obtained by constructing a generating functional, and evaluating the real and imaginary 
parts of the fermion determinant 0]. 
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In this work we follow this second approach, studying the phenomenology of spin 1 
mesons in the framework of the NJL model. We propose a fermionic Lagrangian that shows 
an approximate U{3)l^U{3)r chiral symmetry, explicitly violated by current quark masses, 
as is the case of QCD in the limit of large number of colours Nq- In particular, we revisit the 
predictions given by this model for anomalous radiative decays of vector mesons, taking into 
account the explicit breakdown of the SU{3) flavour symmetry to SU{2) isospin symmetry 
with the assumption = 7^ m^. We perform the bosonization of the NJL theory by 
carrying out an expansion of the fermion determinant in terms of the meson fields. This gives 
rise to a set of one-loop Feynman diagrams from which one can derive the relevant 

effective meson interactions. Taking into account present experimental data, we analyze 
possible input parameters for the model, in order to obtain an acceptable phenomenological 
pattern for the observed anomalous decays of pseudoscalar and vector mesons. 

The paper is organized as follows: In Sect. II we introduce the NJL Lagrangian and the 
bosonization technique. In Sect. Ill we derive the effective Lagrangian to account for vector 
and axial- vector meson interactions, considering the explicit breakdown of flavour SU{3) 
symmetry. This leads to some relations between vector meson masses and decay constants 
that can be written in terms of symmetry breaking parameters. Then, in Sect. IV, we 
concentrate on the anomalous radiative decays of pseudoscalar and vector mesons, which 
proceed through the U{1) axial anomaly. In Sect. V we discuss the numerical results for the 
corresponding branching ratios, in comparison with the present experimental information. 
Finally in Sect. VI we present our conclusions. The Appendix includes a brief description 
of the regularization scheme used throughout our analysis. 



II. NAMBU JONA LASINIO LAGRANGIAN AND BOSONIZATION 



Chiral effective models have proved to reproduce the low energy hadron phenomena 
with significant success. For that reason, a significant effort has been done in order to 
derive these effective schemes from a microscopic theory as QCD. The effective four fermion 
Lagrangian proposed by Nambu and Jona-Lasinio [§, which shows the chiral symmetry of 
QCD, represents one of the most popular ways to achieve this goal. The NJL model provides, 
after proper bosonization, not only the expected effective couplings for scalar, pseudoscalar, 
vector and axial-vector mesons but also the Wess-Zumino-Witten sector which accounts 
for the anomalous meson decays. The explicit form of the NJL Lagrangian reads 



C = q{i ^ ■ 
-2G2 



■mo)q + 2Gi 
1 



(1) 



where q denotes the A^-fiavour quark spinor, 
of the U{N) flavour group (we normalize A° 



0, 



A^^ — 1 are the generators 



'2/N'SL) and stands for the current 
quark mass matrix, which explicitly breaks the chiral symmetry. The coupling constants 
Gi and G2, as well as the quark masses, are introduced as free parameters of the model. 
In the absence of the mass term, the NJL Lagrangian shows at the quantum level the 
SU{N)a ® SU{N)v ® U{l)v symmetry characteristic of massless QCD. 
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It is possible to reduce the fermionic degrees of freedom to bosonic ones by standard 
bosonization techniques. With the introduction of colorless boson fields and an appropriate 
use of the Stratonovich identity, one obtains an effective action where the fermions couple 
to the bosons in a bilinear form. The couplings driven by Gi in (|I]) lead to the introduction 
of scalar and pseudoscalar boson fields, whereas those carrying G2 lead to the inclusion of 
vector and axial-vector mesons. 

To be definite, let us consider the vector meson sector with = 3 flavours, u, d and s. 
By means of the Stratonovich identity, the vector-vector coupling in ([l|) can be transformed 
as 



(2) 



where = —i'^V/^ A"/2. The spin 1 fields can be identified with the usual nonet of 



a=0 



vector mesons, 



V2 Vg 



V 



V2 



72 



V 



K*- 



K*+ 
2ujs 



+ 



^1 



(3) 



which transform in such a way to preserve the chiral symmetry of the original NJL La- 
grangian (and therefore that of QCD). Notice that the first term in the right hand side of 
Eq. (0) is nothing but a mass term for the vector fields V^, thus the vector-meson masses 
are governed by the coupling G2 in the NJL Lagrangian. It can be seen that these masses 
are degenerate in the limit where the quark masses are degenerate. 

Now, the quark fields can be integrated out, leading to an effective Lagrangian which 
only contains bosonic degrees of freedom. This procedure can be carried out by taking into 
account the generating functional 



Z = N j VVVqVq exp |i j dx 



4G. 



-Tr + q{i ^ — rho + i y)q 



(4) 



and performing the calculation of the fermion determinant (a detailed analysis can be found 
in |J7[). A similar procedure can be followed for the full NJL Lagrangian (p, leading to the 
interactions involving scalar, pseudoscalar and axial-vector bosons. In this way, the final 
effective Lagrangian is written only in terms of spin and spin 1 colorless hadron fields. 

Here we have performed the bosonization by carrying out an expansion of the fermion 
determinant, which gives rise to a set of one-loop Feynman diagrams [§. In order to build 
up the final effective meson Lagrangian, we have obtained the local part of the relevant 
interactions by taking the leading order of a gradient expansion in powers of the external 
momenta, and considering the dominant contributions in 1/Nc- The divergent integrals 
have been treated using a proper-time regularization scheme with a momentum cut-off, 
keeping the leading contributions. The procedure is developed in detail in the next section. 
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III. EFFECTIVE MESON COUPLINGS AND SYMMETRY-BREAKING 

PARAMETERS 



In this section we describe the effective meson interactions derived from the NJL La- 
grangian. We begin by calculating the relevant one-loop self-energy diagrams contributing 
to the meson kinetic and mass terms, and those leading to the effective weak and strong 
meson decay couplings. Notice that the meson fields introduced through the bosonization 
technique acquire their dynamics via the self-energy one-loop diagrams. 

In our framework, the poles of the quark propagators occur at the so-called constituent 
masses, i.e. djTiamically generated masses which arise as a consequence of the spontaneous 
breakdown of the chiral symmetry. The constituent masses appear as solutions of Schwinger- 
Dyson self-consistency equations (gap equations), which are governed by the scalar coupling 
Gi in the NJL Lagrangian (|I]). Here we will not concentrate in the scalar sector of the 
model, thus the constituent quark masses, as well as the pseudoscalar boson masses, will be 
considered as free parameters. 

Our main interest is focused in the anomalous radiative decays of vector mesons. Never- 
theless, we need to consider also the axial-vector sector of the model, since the axial-vectors 
mix with the pseudoscalars at the one-loop level [P — A mixing). 



A. Kinetic terms and masses of the spin-1 mesons 



The generating functional (|) gives rise to effective kinetic terms for the spin-1 vector 
mesons via one-loop diagrams, as shown in Fig. 1(a). The analysis in the case of the axial- 
vector mesons can be performed in an analogous way, with the additional ingredient of P — A 
mixing given by Fig. 1(b). 

The loop in Fig. 1(a) gives a contribution to the vector meson self-energy given by 



iN, 



c 



rTr 



^— mi 



(27r)4 (k-py-ml '^k'^-ml 



(5) 



where mi and m2 are the constituent masses of the quarks entering the loop, A''^; is the 
number of colors, and the trace acts over the flavour and Dirac indices. 

As stated above, we will take only the leading order in the external momentum p, which 
means to evaluate the integral at p = after extracting the relevant kinematical factors. In 
this case, this is equivalent to consider only the divergent piece of (^: 



nj.^^ = h{mi,m2) 



mi. 



where 



l2{mi,mj) = 



—I 



c 



(27r) 



(A;2 — mf){k'^ — m^) 



(6) 



(7) 



As it was previously mentioned, in order to regularize the divergence we use the proper-time 
regularization scheme [jl3[ with a cut-off A, which will be treated as a free parameter of the 
model. Details of the procedure can be found in the Appendix. We obtain 
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From (P), the kinetic terms for the vector mesons in the effective Lagrangian are given 

by 



1 2 

- - J2(,m„, m„j [p^^p^" + 2p^^, 
+2/? (k*I^K*'^'' + , (9) 

where V"'^^ = 9'^^'^ - S^V'^, and 

_ l2{ms,ms) _ hjmu^ms) . . 

l2[mu,mu) l2{mu,mu) 

parameterize the magnitude of the SU (3) flavour symmetry breaking. 

The kinetic Lagrangian in (|^) has been expressed in terms of the vector fields in @, 
with the additional rotation 

tUg = cos 9o + uj sin 6^0 

cui = — sin 6*0 + L<j cos 6*0 , (11) 

which diagonalizes the neutral sector. It is easy to see that here the rotation is "ideal", i.e., 
the spin 1 mass eigenstates p and uj are composed by pure light u and d quarks, while the 
(j) meson is a bound state ss. The a; — rotation angle is given by sin6'o = Slight 
deviations from the ideal mixing condition will be considered below to allow for the decay 
(f) —>■ TT^'-f observed recently. 

The mass terms for the vector mesons are given by the (V^)^ term in (^, plus a divergent 
one-loop contribution given by the second term in the square brackets in (j^), which vanishes 
in the SU{3) flavour limit. This leads to 

+ {^,-m^fPh{m^,m^){K*lK*~^ + K*lK*'^) . (12) 

Notice that (as expected) the mass terms turn out to be diagonal in the (cu, 0) basis. 

We proceed now to the wave function renormalization required by the kinetic terms in 

1 /2 

(y). The vector meson fields can be properly redefined by — > Zy V^, with 

2 

= Z-^ = - hirriu, rriu) (13a) 
Z^l = '^P h{mu, m„) = P Z;' (13b) 
Z^'^ ='^al2imu,mu) = aZ~^ . (13c) 



Then from (O) one obtains 
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m. 



4G2 



3 



(14) 



thus the meson mass can be written in terms of the p mass and the flavour symmetry 
breaking parameter a. In the case of the K* , the corresponding mass relation includes both 
the parameter (3 and a quark-mass dependent contribution that arises from the loop in Fig. 
1(a). 

A similar analysis can be performed for the axial-vector meson sector. By replacing 
7a — ^ 7a 75 in the quark-meson vertices in Fig. 1 (a), one finds 



leading to the kinetic Lagrangian 



3 (m2 + nil] 



(15) 



{A) 
kin 



1 2 



(16) 



and the mass terms 
1 



a.^i + 2a+ a^^ + /i^/f + /' J'^ + 2Kt K^'' + 2KIk[^ 



+ (m. + m,)^/3(iri+iri-^ + KilK'/) + 2miaf[J"t 



(17) 



As in the vector meson case, the Lagrangian is diagonal in the chosen basis, where /i and 
f'l are obtained from the U{3) states Ag and Ai through an ideal rotation. Notice that 
both neutral and charged axial-vector boson masses receive a contribution proportional 
to G^^, which arises from field transformations as in Eq. (^, plus additional (positive) 
contributions proportional to the quark masses. In this way, the axial-vector mesons are in 
general expected to be heavier than their vector meson counterparts. 

Finally, we have to take into account the mixing between the axial-vector mesons and 
the pseudoscalars. This requires the analysis of the pseudoscalar-axial vector couplings 
and the pseudoscalar kinetic terms, given by the one-loop diagrams in Fig. 1 (b) and (c) 
respectively. The effective couplings read 



^mix = -ihifUu, m„) [2m„(9^7r°< + S^vr+ai ^ + d^n a^^ + d^'quf ) 

+ 2amsd^7^J'^] , 
while the pseudoscalar kinetic Lagrangian is given by 



(18) 



{P) 

kin 



1 



+ 2pdf,K+d''K- + 2pd^K°df'K^ 



(19) 
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Once again, we have chosen a basis for the neutral fields in which the states rju, rjs are 
obtained from the U{3) states t]s, rji through an ideal rotation. However, these states cannot 
be treated as approximate mass eigenstates due to the presence of the U (1)^ anomaly, which 
breaks the ?7(3) symmetry down to SU{?)). Although formally of order 1/Nc, this anomaly 
leads to a relatively large mass splitting between the observed r] and rj' physical states. 
Still, large Nc considerations are shown to be powerful enough to deal with the interactions 
of the rji and rjs fields as members of an U{3) nonet. As stated, we will not concentrate 
here in the pseudoscalar mass sector, and consequently, in the U{1)a symmetry breaking 
mechanism responsible for the rj' mass. In the NJL framework, one way to proceed is to 



include the so-called 't Hooft interaction |I^p!5[ , which emulates the effect of the anomaly. 



Instead of following this way (which means to include six-fermion interaction vertices) we 
will introduce the rj — rj' mixing angle as a parameter of the model. Thus we write rju and 
rjs in terms of the physical states as 

rju = cos ifp Tj + snufp rj' 

rjs = — simp p rj + cos Lpp rj' . (20) 

In order to write the effective Lagrangian in terms of the physical fields, one needs not 
only the relevant wave function renormalizations but also the diagonalization of the P — A 
couplings. This can be achieved by means of the transformations 

P^Z]I'P, A^^Z'/A^ + CpZ]I'd^P, (21) 

where P = n, K, rju, rjg, and 

= i2l2{mu,mu)Zp^^ = C 
Ck = i2 hirriu, rriujZp — - — ^ 

771 

Cr,^ = i 2 hiniu, mu)Zp^ . (22) 

rrij, 

The wave function renormalization for the pseudoscalar sector yields 



Z-.' = Z;,: = h{mu,mu)[l-'^-^ 



Zk = Phimu^mu) 



2 ml. 



Z^/ = ahiiTLu, rriu) ( 1 - ) , (23) 



where the axial- vector meson masses can be read from (|T^). As can be immediately seen 
from (|16|) , the wave function renormalization for the axial-vector mesons is equivalent to 
that performed in the vector meson sector, c.f. Eq. (|13|). Using these relations and those in 
(|l^), the axial-vector meson masses turn out to be constrained by 
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2 2 I p 2 

= m|-. + QmuTUs (24) 
m J, = m| + 6ml . 



B. Weak and strong decay constants 



We want to analyze the ability of the effective model under consideration to describe the 
low energy behaviour of pseudoscalar and vector mesons. Let us begin by studying some 
fundamental decay processes, in order to evaluate the relevant decay constants in terms of 
the model parameters. 

The pseudoscalar weak decay constants and Fk can be obtained from the one-loop 
transitions represented in Fig. 2, where P = tt, K. As in the self-energy case, these diagrams 
are found to be logarithmically divergent and can be regularized using the proper-time 
scheme. We will keep as before only the leading contributions in the limit of vanishing 
external momenta, which amounts to take only the divergent piece of the loop integrals. 

The decay constants Fp are defined as usual by 



(0|J^5(0)|P(g))=«v^Fpg^. 

It is easy to see that the contributions from Fig. 2 (a) and (b) lead to 

F^ = 2 hiruu, mu)m^Zl/'^ (1 + z 2 m„C^) 
= 2/2(m„,mJm,Zy2 _ 

F, = iMu. m.) /3 (m„ + m,) [1 + z (m„ + m,)C,] 

3 (m„ + ms)2 



(25) 



h{mu, mu)l3 {ruu + ms)Z 



1/2 
K 



1 - 



(26) 



where have made use of Eqs. (|2^) . In these equations, the first terms arise from the graph 
in Fig. 2 (a), whereas the contributions proportional to Ct^^k are due to the P — A mixing 
diagrams. Using now Eqs. (|23|) , we end up with the simple relations 



2m„ 
7I/2 

Zj-k 



(27) 



Next, let us study the basic decay constants Qp, qk* and g^, which account for the strong 
decays of the vector mesons into two pseudoscalars. The decay constants can be defined 
from the total decay rates according to 



TTTT 



r(p 

T{K* Kn) 



9l_ rrip 
An 12 



9k- ^k* 



4m: 



3/2 



47r 16 



{rriK + m^) 



1 - 



3/2 
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KK) = r(0 ^ K^K~ 



5^ 

A-K 24 




m3. 



Notice that for the — > KK width we have taken into account the mass difference between 
the neutral and charged kaons. Though tiny, this mass difference becomes important in view 
of the very narrow phase space allowed. 

Within the effective model analyzed here, the decay rates (^) are dominated by the 
one-loop graphs shown in Fig. 3. As before, we keep here the dominant, logarithmically 
divergent contribution, and the leading order in the external momenta (the corresponding 
diagrams with two axial-vector legs vanish in this limit). From the evaluation of the p ^ tttt 
and p — > aiTc diagrams, and using (|^), (0) and (^3]), we find 



g, = him^, m^)Zy^Z^ (1 + z 2 m„ a) = . (29) 

In the case of the evaluation of the dominant contribution is more involved, since the 
loop includes two different quark propagators. The loop integral has the form 



. Nc 

-I 



(27r)4 J ^ ^ (P - ml) (A;2 - mlf ' 

Once again the divergence is logarithmic, thus it is natural to express the integral in terms 
of himi, rrij). However, owing to the flavour symmetry breakdown, there is an ambiguity at 
the moment of choosing the infinite part of (PD]). Here we follow the criterion of respecting 



the natural hierarchy given by the quark content of the decaying vector meson, that means, 
we use the k"^ factor in the integrand to remove the pole. We have then 

Is a/2(m„,m„) . (31) 

Taking into account the wave function renormalization for the and K mesons, together 
with the corresponding pseudoscalar-axial vector mixing, we end up with 



a 

p- 9p ■ (32) 

In the same way we proceed for the case of the K* to two meson decay. Here, due to the 
SU (3) symmetry breaking, the analogous of Eq. ( p9D for g^* leads to a relation between gx* 
and gp which can be written in terms of Z^^ k and Ct^^k- At first approximation, however, 
the dependence on these parameters cancels and one can write 

Qk' ^ gp ■ (33) 

In this way, both (p and K* strong decay constants can be predicted in terms of the 
measured value of gp and the symmetry breaking parameters a and (3. The corresponding 
numerical evaluation will be done in Sect. V. 
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IV. ANOMALOUS RADIATIVE DECAYS OF PSEUDOSCALAR AND VECTOR 

MESONS 



Our goal is to test the ability of this relatively simple, NJL-like model to reproduce the 
observed pattern for the anomalous radiative decays of pseudoscalar and vector mesons once 
the breakdown of chiral symmetry is taken into account. These decays are mainly originated 
by the U{1)a anomaly. In the framework of our model, they occur through the well-known 
triangle graphs shown in Fig. 4, which lead to an anomalous Wess-Zumino-Witten effective 
Lagrangian. The corresponding diagrams including an axial-vector which mixes with the 
pseudoscalar meson gives a vanishing contribution This can be seen e.g. by looking at 
the effective anomalous Lagrangian obtained by gauging the Wess-Zumino-Witten action, 
as shown in Ref. ||^. 

The couplings between the vector mesons and the outgoing (on-shell) photons follow 
from the assumption of vector meson dominance of the electromagnetic interactions. In our 
framework, these couplings can be obtained from quark-loop diagrams connecting photons 
and vector mesons (see Ref. [0,0 for details). After the proper inclusion of the symmetry 
breaking parameters a and (3, the V — 7 couplings read 



rr = - < z;^" [p. + \-. + ^^v^ ) • (34) 



p / 

Let us begin by analyzing the decays of the neutral pseudoscalars 7r°, 7] and rj' into two 
photons. These processes arise at the one-loop level from the triangle diagrams shown in 
Fig. 4 (a). The diagrams are found to be convergent, and the limit of vanishing external 
momenta can be taken trivially. For vr*^ 77 we get 



r(-"-77) = ^^5^j (35) 

where a is the electromagnetic fine structure constant. This expression coincides with the 
result obtained in standard Chiral Perturbation Theory [|l| . In the case of rj and rj' decays to 
two photons, one has the additional problem oi rj — t]' mixing. In our framework, instead of 
dealing with the f/(3) states r/i and rjs, it is natural to work with the rju and rjs states defined 
above, which determine the flavour content of the quark propagators in the loop. Taking 
into account the corresponding pseudoscalar wave-function renormalizations, we find 



647r3 



Y(7]' 77) - " 



647r3 



V2 . 

cos 09 p Sm 09 p 

■sm(pp + — — cos (p p 



(36) 



where, in analogy with and Fk, we have defined 



^^^^2- (37) 
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Now we analyze the radiative vector meson decays V 
can be conveniently parameterized as 



P7. The corresponding widths 



T(y P7) 



'^9p Cyp 

24^ ~fT 



My 



(38) 



where the coefficients Cyp can be obtained from the loop diagrams shown in Fig. 4 (b). 
Clearly, the evaluation of the relevant triangle quark loops is entirely similar to that per- 
formed for the P — >■ 77 case. 



For the processes p°, 



UJ 



7 we can use Eqs. (|13a]) , (^) and (|29| ) to eliminate the 



pseudoscalar and spin 1 wave function renormalization factors. In this way the coefficients 
CpO^o and C^T^o are simply given by 



CpO^O = 1 



a 



(39) 



The same procedure applies in the case ofV rj, processes, with the additional inclusion 
of the 1] — rj' mixing angle and the flavour symmetry breaking parameter a. We have 



3 cosipp 

COSipp 

2P^sin ipp/{y/aFs] 
2P^cos Lpp/{^/aFs 



(40) 



where ipp is the rju — f]s mixing angle defined in (pOD . 

For the decays K* K'j, the calculation is shghtly more complicated due to the presence 
of different constituent masses for the quark propagators in the loop. The coefficients depend 
now on the symmetry breaking parameter j3 and the mass ratio A = rris/mu according to 



Ck*+k+ 



f{X)Fj{^FK) 
g{X)Fj{^FK) 



(41) 



where functions /(A) and g{\) are defined as 

■I + 6A ^2 



/(A) 



2(1- A2) 
A In A^ 



A^ ^ A(2A2 + l)lnA^ 



;i-A2)2 



1 + 



(42) 



Finally let us consider the decays rj' — > p7 and r]' — > 077, which can be treated in a 
completely similar way just taking into account the spin character of the decaying particle 
when averaging over the initial spin states. Using the same parameterization as in (^) we 
obtain 



CpOrj' 
Cuirj' 



3 sin 

sinifip 



(43) 



We have skipped in this discussion the decay 7r°7 since the corresponding coefficient 
C<A^ vanishes in our framework, owing to the purely strange flavour content of the meson. 
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However, in view of the present experimental results, it is interesting to remove the ideal 
mixing angle condition in the u — (f) sector allowing for a tiny nonstrange component for the 
field. If this is parameterized by a small mixing angle (fv, one trivially gets 



C^T,o = C^^o sin v^y = 3 sin (pv ■ (44) 

The present experimental value for the branching ratio — * tt^'j leads to a mixing angle ipv 
of about 3.2° [|T^,|r^, which does not represent a significant change in the other, nonvanishing 
radiative decay widths involving and uj. 



V. MODEL PARAMETERS AND NUMERICAL ANALYSIS 

In order to perform the phenomenological analysis of the model, and to obtain definite 
predictions for pseudoscalar and vector meson decay processes, it is necessary to establish a 
strategy, choosing a suitable set of parameters to be used as input values. 

Let us begin with the basic decay processes analyzed in Section III. Using the vr^ and K"^ 
weak decays, it is possible to write the pseudoscalar wave function renormalization factors 

and Zk in terms of the well-measured decay constants F^^ and F^, which will be taken 
as input parameters. In the same way, the wave function renormalization factor Zp can be 
obtained from the phenomenological value ~ 6 arising from (|28|) . Now, from (H) and 
( |13a|) , we can relate gp with the loop integral J2(m„,m„) to estimate 

^ ~ 0.26 . (45) 



Next we use the mass relations (14) to determine a phenomenological acceptable set of values 



for the cut-off A and the strange quark mass m^. A reasonable choice is 

rris ^ 510 MeV , A ~ 1.2 GeV , (46) 

which lead to = 0.99 GeV (exp. 1.02 GeV) and ttik* = 0.91 GeV (exp. 0.89 GeV), with 
a light quark mass ruu = rrid — 310 MeV. The symmetry breaking parameters a, j3 and A 
defined in the previous sections are then given by 

a ~ 0.6 , /? ~ 0.76 , A ~ 1.6 . (47) 

Using the relations (|3^) and (|33D , the phenomenological value of Qp, together with the 
results in ( ^7D allow to predict the values for the strong decay constants qk* and g^. In this 
way we obtain 

T{K* Kit) ^ 44 MeV 

r(0 ^ K^K°) ~ 1.3 MeV (48) 
r(0 ^ K+K') ~ 2.0 MeV 

while the corresponding experimental values read r{K* — > Kit) = 50.8 ± .9 MeV, r(0 
K^K^) = 1.51 ± 0.03 MeV and rf^ ^ K+R-) = 2.19 ± 0.03 MeV HI. The results are in 



agreement with the predicted values within an accuracy of a 15% level, which can be taken 
as a lower bound for the expected intrinsic theoretical error of the model. 
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As an alternative procedure, we could have derived the values for F^^ and Fk from the 



relations (23), where the axial- vector meson masses can be taken either from Eqs. (24) or just 



as input parameters. For instance, with the chosen value of m„ we can obtain from Eq. ( p^ 
a prediction for the ai mass of about 1080 MeV [0. This value is somewhat low, but once 



again agrees within a 15% accuracy with the present experimental result of = 1230 ±40 
MeV. In fact, it is natural to expect the model to describe the axial- vector meson sector only 
in a roughly approximate way, due to the relatively high masses involved, the broad character 
of the resonances and the admixtures with close states with same quantum numbers. Hence 
we choose to rely on the well-known low energy decay constants and Fk as our input 
magnitudes. 

Finally, using the analytical expressions obtained in the previous section, we can look at 
the model predictions for pseudoscalar and vector meson radiative decays. In the case of the 
decays involving t] and r]' mesons, it is necessary to fix two more input parameters, which 
can be chosen as the decay constant Fg (or equivalently the pseudoscalar wave function 
renormalization Z^^) and the rji — rjg mixing angle 6 p. The latter is related with v?p by 

V^p = ^ - ^0 + , (49) 

where 6*0 — 35.3° is the "ideal" rotation angle introduced above. 

Our numerical results are shown in Table I, where we quote the values corresponding 
to four different sets for Fg and 6p, together with the present experimental widths [1^. In 



order to get the theoretical results, we have used the physical values for the pseudoscalar and 
vector meson masses entering Eq. (|38|). We have chosen the values of Fs and 6p that lead 
to a better agreement for the whole set of decays, always having in mind that an intrinsic 
theoretical error of at least 15% can be expected from the NJL scheme. In general, we find 
that the best results for 6p lie within the usually expected range between —10° and —20° 
|T^,|T9|, while the model seems to prefer a rather large value for the ratio Fs/F^^, higher than 
3/2. We quote for completeness the results corresponding to a ratio Fs/ F^, = 1.4, favored by 
large- A''^; arguments , which leads in this case to too large widths for the decays involving 
7] and mesons. 

When varying Fg and 6p, and comparing with the experimental values, it is found that an 
accurate agreement for certain widths worsens the predictions in other cases. Nevertheless, 
considering the relatively important intrinsic theoretical error of the model, it can be said 
that the values in Table I show a phenomenologically acceptable pattern for the pseudoscalar 
and vector meson anomalous radiative decays. Our results are found to be of similar quality 
to those presented e.g. in Ref. where the authors start from a nonlinear sigma model, 
and Ref. |jl2| , where a different parameterization of the chiral-fiavour symmetry breaking is 
introduced (notice that several experimental values have changed in the last years). 



VI. CONCLUSIONS 

Starting from an effective four-fermion NJL Lagrangian, we have analyzed the phe- 
nomenology associated with anomalous radiative decays of pseudoscalar and vector mesons. 
The effective interactions between spin and spin 1 mesons have been derived from the 
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bosonization of the fermionic theory through the evaluation of the relevant quark loop dia- 
grams, taking the leading order both in the external momenta and I/Nq power expansions. 
The divergent loops have been treated using a proper-time regularization scheme with a 
momentum cut-off A. 

We have taken into account the breakdown of the chiral symmetry considering constituent 
quark masses m„ = 7^ irts. The departure of the effective couplings from the global 
SU{2>) flavour symmetry limit has been introduced through the parameters a and /3, which 
correspond to properly regularized fermion loop integrals. 

In order to perform the phenomenological analysis, we have chosen a few fundamental 
input parameters, namely the basic decay constants F^r, and Qp, the p meson mass, the 
momentum cut-off (or equivalently the light quark mass m^), and the constituent mass of the 
strange quark nis- We have also taken into account the physical values of the pseudoscalar 
meson masses. To account for the decays involving rj and rj' mesons, this set has been 
complemented with the parameter Fg and the 7] — 7]' mixing angle 6p. We have shown that 
with the choice A ~ 1.2 GeV and ~ 510 MeV this simple model leads to an acceptable 
pattern for the vector meson mass spectrum, as well as the branching ratios for the strong 
decays p —>■ tttt, cf) —>■ KK and K* —>■ Kn. The agreement is found within a 15% level, that 
can be taken as a lower bound for the intrinsic theoretical error of the model. 

We have evaluated within this framework the branching ratios for anomalous radiative 
decays of pseudoscalar and vector mesons. In the case of those processes involving the 77 
and r]', the contrast between the predictions of our model and the present experimental 
data is optimized for a 6p mixing angle between —10° and —20°, in agreement with usual 
expectations. On the other hand, we find that the preferred ratio Fs/Ft^ lies in the range 
1.75 to 2, which is a rather high value in comparison with the result arising from large Nc 
considerations. In general, taking into account the relatively large theoretical error, we see 
that the model leads to reasonably good predictions for the main radiative pseudoscalar and 
vector meson decay widths. The quality of our results is comparable to that obtained in 
previous works which use different effective models and/or parameterizations for the flavour 
symmetry breaking effects. 



The divergent integrals /2(^i, rnj) have been regularized within the proper-time scheme 



which holds for A sufficiently large. If we let A be a momentum cut-off for our theory, we 
have for equal masses = rrij 



APPENDIX: 



P^ . One makes use of the relation 




(Al) 




(2vr)4 



Nc 



I 



1/A2 



00 



ds s e "'^^B 



IGtt^ 



Nc 




) 



(A2) 
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where fc^ is the momentum in Euchdean space, and r{a,x) is the Incomplete Gamma 
Function, 

T{a,x) = / e"* dt . (A3) 

Jx 



In the case where rrii ^ rrij, Eq. (^) can be generalized using the Feynman parameter- 
ization 

1 1 

dx —2 „, , (A4) 



(4 + m^){kl + m]) Jo {kl + B\x)y ' 

where B'^{x) = mf + (m^ — mf) x. Using the same regularization procedure as before one 
ends up with 

rfxr[o,^) . (A5) 
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(a) 



(b) 



(c) 



FIG. 1. (a) One-loop self-energy diagrams for vector and axial vector mesons, (b) pseu- 
doscalar-axial vector mixing diagram, (c) pseudoscalar self-energy diagram. 




(a) (b) 
FIG. 2. One-loop diagrams accounting for weak decays of pseudoscalar mesons. 





(a) 



(b) 



FIG. 3. One-loop diagrams accounting for strong decays of vector mesons into two pseu- 
doscalar s. 
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(a) 



(b) 



FIG. 4. One-loop diagrams accounting for anomalous decays of (a) pseudoscalars and (b) vector 
mesons 
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TABLES 



TABLE I. Radiative decays of pseudoscalar and vector mesons. The second and third columns 
show the model predictions for different values of the ratio Fg/F-,^ and the r/i — ryg mixing angle Op. 
The present experimental values for the decay rates are quoted in the last column. 



Process Vnjl (KeV) Te^p (KeV) 

Fs/F^ = 2 Fs/F^ = 1.75 = 1.75 F,/F, = 1.4 





= -10° 


0p = -15° 


Op = -20° 


0p = -20° 




7r° — > 77 


7.7 X 10-^ 


7.7 X 10-3 


7.7 X 10-3 


7.7 X 10-3 


(7.7 ±0.6) X 10 


77^77 


0.53 


0.63 


0.76 


0.71 


0.46 ± 0.04 


T]' 77 


4.97 


4.48 


3.79 


4.16 


4.28 ± 0.44 


rj' p'^7 


108 


89.5 


71.1 


71.1 


59 ±5 


ij' — > co^ 


10.0 


8.2 


6.5 


6.5 


6.1 ±0.8 




85 


85 


85 


85 


102 ± 25 


o; — > 7r^7 


806 


806 


806 


806 


720 ± 40 


P° ^ '77 


51 


60 


68 


68 


36± 14 


(J ^ 777 


6.6 


7.7 


8.8 


8.8 


5.5 ±0.8 


^ ^ ?77 


65.0 


70.0 


55.6 


86.9 


57.8 ± 1.5 


0-^ ?7'7 


0.30 


0.46 


0.52 


0.82 


0.30 ±0.16 




62 


62 


62 


62 


50 ±5 




160 


160 


160 


160 


117 ± 10 
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